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The exotic nature of many strongly correlated materials at reasonably high temperatures, for instance cuprate 
superconductors in their normal state, has lead to the suggestion that such behavior occurs within a quantum 
critical region where the physics is controlled by the influence of a phase transition down at zero temperature. 
Such a scenario can be thought of as a bottom-up approach, with the zero temperature mechanisms finding a way 
to manifest critical behavior at high temperatures. Here we propose an alternative, top-down, mechanism by 
which strong kinematic constraints that can only be broken at extremely high temperatures are responsible for 
critical behavior at intermediate but still high temperatures. This critical behavior may extend all the way down 
to zero temperature, but this outcome is not one of necessity, and the system may order at low temperatures. We 
provide explicit examples of such high-temperature criticality when additional strong interactions are introduced 
in quantum Heisenberg, transverse field Ising, and some bosonic lattice models. 



I. INTRODUCTION 

Strongly correlated systems display very rich phase dia- 
grams upon varying thermodynamic parameters such as tem- 
perature, pressure, doping concentration, magnetic field, and 
so on. 1 At sufficiently low temperatures, distinct phases of 
matter appear that are characterized by (sometimes coexist- 
ing) long-range order of, say, the magnetic, charge, orbital, 
or superconducting type. When transitions between different 
zero temperature quantum phases are continuous, fluctuations 
occur at all length scales and lead to power law behavior for 
correlation functions of the order parameter at the critical cou- 
pling. Tuning the temperature slightly away from zero de- 
creases the strength of the order parameter fluctuations but as 
long as they remain sufficiently strong they lead to scaling 
laws that are insensitive to the microscopic details and, to a 
large degree, universality has emerged. 

The temperature range for which scaling laws apply is a 
measure of how strong fluctuations are. In strongly correlated 
systems such as organic materials, high-T c superconductors, 
etc, these scaling laws extend to surprisingly high tempera- 
tures above the critical temperature. This fact is attributed to 
the presence of very strong one- or two-dimensional fluctua- 
tions of the order parameter that are predominantly quantum 
at sufficiently short length scales. How large the tempera- 
ture is at which quantum fluctuations are effectively observed 
is a matter of intense debate. 2 For example it has been pro- 
posed that many exotic properties in the so-called pseudogap 
regime of high-T c superconductors originate in a hidden quan- 
tum critical pointi^i^ In this scenario, increasing the doping 
concentration induces at zero temperature a continuous phase 
transition between two different states of matter below the 
superconducting dome, that is reflected in the strange metal- 
lic properties of the normal state above the superconducting 
dome. 

In this paper, we present an alternative picture where 
the critical behavior in some strongly interacting systems is 
present at high temperatures due to strong kinematic con- 
straints in the quantum Hamiltonian. The system remains 
critical as long as the constraints are respected, i.e., below 



some large energy scale, corresponding to the largest coupling 
in the Hamiltonian. This critical phase becomes the physi- 
cally important, universal feature from which a scaling regime 
can descend. The zero-temperature physics is instead system- 
specific, with a rich variety of different ordered phases, and 
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FIG. 1 : (Color Online) - Generic phase diagram for a strongly con- 
strained quantum system that satisfies conditions 1) to 6) in Sec.Q 
The parameter space encodes the competition between two energy 
scales, the temperature T and a characteristic energy scale g that se- 
lects a classical ordered phase, while F eff and U are held fixed. A 
quantum critical scaling regime, if it exists, is restricted to a region 
represented by the upper half of a disk of radius ~ |r cff j and cen- 
tered at the origin (g, T) — (0, 0) of parameter space. The focus of 
this paper is on the constrained entropic scaling regime at tempera- 
tures intermediate between the small characteristic quantum energy 
scale r cff and the large characteristic energy scale U set by a strong 
constraint. At a fixed temperature, the constrained entropic scaling 
regime terminates in a phase transition that needs not be continu- 
ous upon increasing \g\. At fixed g, the constrained entropic scaling 
regime crosses over to the conventional high-temperature phase, say 
a paramagnetic one for spin degrees of freedom, when T is of the or- 
der of U. The transition from the constrained entropic scaling regime 
to the quantum regime upon lowering T at fixed g is system specific. 
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may or may not allow for the finite temperature criticality to 
survive all the way to zero-temperature. We also discuss how 
this picture brings about four distinct scenarios for the finite 
temperature behavior of the system, depending on the type of 
interactions present in the quantum Hamiltonian. In particu- 
lar, we show how one scenario naturally leads to a quantum 
system exhibiting an exotic correlation length that increases 
with increasing temperature, over a wide range of tempera- 
tures. 

The making of a system displaying high-temperature criti- 
cality is as follows. 

1) At least one characteristic energy scale in the quantum 
Hamiltonian, say U, is much larger than all others. Thus, U 
splits the Hilbert space in sectors separated by large energy 
gaps. In particular, the infinite U limit projects the Hilbert 
space onto the space Tl v of allowed states. 

2) The quantum dynamics is generated by the terms that 
do not commute with the [/-term in the quantum Hamilto- 
nian. However, one needs a process of order n > 2 in the 
perturbative expansion in the characteristic energy scale of 
these terms, V, to generate an effective quantum coupling be- 
tween allowed states. Hence, this effective quantum coupling 
r cff ~ TiT/UY 1 - 1 can be very small, |r off | < |T|. 

3) By squashing down the effective quantum coupling T eff , 
one opens a hierarchy of temperature scales, which we dis- 
cuss in detail in the paper. In particular, already at reasonably 
low temperatures (above the small |r cff |), the effective quan- 
tum term can be neglected in the calculation of equilibrium 
thermodynamic quantities. 

4) Although the equilibrium thermodynamics for T ^s> \T cS \ 
is classical, the dynamics is still quantum. The reason is the 
following. If the system were to rely on thermally activated 
processes to move within the restricted Hilbert space H v , 
the characteristic time-scale would be, for a small system- 
bath coupling 7 > 0, t t = 7 _1 exp(U/T), which is astro- 
nomical for temperatures well below the large energy scale U. 
There are also virtual processes due to system-bath coupling 
that bypass thermal activation, but whose effective coupling is 
suppressed: 7 off ~ 7(7 /U) n ~ 1 <C 7, much alike the intrin- 
sic terms that do not commute with the [/-term. The waiting 
time for a quantum tunneling event is Tq = min(|r^ | , 7^), 
which depends only algebraically on U. Since 7 <C |T| (the 
bath coupling should be the weakest term in order not to per- 
turb the energy levels of the system), the intrinsic quantum 
dynamics provides the smallest dynamical time scale. Hence, 
a phantom of quantum mechanics in the form of sporadic tun- 
neling events between which coherence is lost provides the 
fastest mechanism for the system to reach classical thermody- 
namic equilibrium when \T cS \ «T <C U. 

5) The strong constraint on the allowed states imposed by tak- 
ing the limit U — > 00 first makes the system critical in the 
constrained entropic limit T — » 00 with all the remaining 
characteristic energy scales held fixed. We call such a criti- 



cal point a constrained entropic critical point. The system's 
properties at low temperatures compared to U are controlled 
by the close proximity to this purely constrained entropic crit- 
ical point. Ice, coloring, and dimer models provide examples 
of constrained entropic critical points in systems with hard 
constraints. 

6) Finally, if there are other terms in the Hamiltonian with 
characteristic energy scale g, \g\ <C U, but that commute 
with the U -term, there can be different phases and transitions 
among them when T/\g\ is 0(1). When U > T > \g\, 
the physics is controlled by the proximity to the constrained 
entropic critical point. Even though the constrained entropic 
critical point may be unstable, the nearby renormalization 
group (RG) trajectories feel its presence until the RG scale 
T ~ U is reached, beyond which the constraint becomes im- 
material. The featureless unconstrained (stable) fixed point 
takes over at that stage. 

The phase diagram in Fig. \l\ exemplifies the high- 
temperature critical behavior arising from kinematic con- 
straints that we highlight in this paper. It shows a large re- 
gion - the constrained entropic critical regime - standing in 
between two ordered phases. The constrained entropic criti- 
cal region exists because of the strong constraints imposed by 
the large energy scale U in the problem, and it covers a high- 
temperature range that ends at the extreme limit of tempera- 
tures of order U, where the system becomes featureless in that 
it is controlled by its proximity to the unconstrained entropic 
fixed point (paramagnetic phase). The ordered phases to the 
left or right of the constrained entropic critical regime have a 
classical origin in Fig.^ since there the second largest energy 
scale \g\ is associated to an operator that commutes with the 
constraint. If we fix r off and vary g, as in Fig. ^ the high- 
temperature constrained entropic critical region sits on top of 
a quantum phase with radius of the order of |T eff | around the 
origin T = g = 0. Quantum criticality may exist only within 
this small region. 

We illustrate this constraint-based critical behavior at high 
temperatures using simple case studies, a constrained quan- 
tum Heisenberg model and a constrained Ising model in a 
transverse field both on the honeycomb lattice, and a con- 
strained bosonic model on the square lattice that leads to the 
quantum dimer model at low energies. The single-band Hub- 
bard model, where a large onsite repulsive term U is the dom- 
inant energy scale, fails to fulfill condition 5), and correlations 
decay exponentially fast with separation beyond a character- 
istic length scale of the order of the lattice spacing. However, 
this situation may change if one considers extended Hubbard 
models with large nearest-neighbor coupling V, next nearest- 
neighbor coupling V' , etc, at commensurate fillings. We shall 
comment on this situation in the conclusions, and discuss, for 
instance, the possible connection between the ideas of con- 
strained entropic criticality and those of fluctuating stripes. 5 
In addition, power-law behavior at high temperatures was also 
shown to occur in the context of locally fluctuating bond cur- 
rents in d-density wave states by Chakravarty in Ref. |y, who 
also recognized the importance of the constraints in determin- 
ing the long distance correlations in the system. 
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The plan of the paper is as follows. Four examples of 
constrained quantum models are introduced in Sec. |H] The 
regimes of temperature that are of relevance to this paper 
are presented in Sec. |in| The constrained entropic scaling 
regime is described in Sec. IIVI Realizations of constrained 
entropic scaling regimes for quantum XXZ Heisenberg and 
transverse field Ising models on the one hand, and lattice 
bosonic models on the other hand, can be found in Sec. fvl 
and Sec. lVIl respectively. In Sec. lVIII we discuss the regime 
when the quantum scale |T eff | is the largest scale below U, 
i.e., \g\ < |T cff | < U, and we conclude in Sec. lVTTTl 



II. FOUR EXAMPLES OF CONSTRAINED QUANTUM 
MODELS 

Constrained quantum mechanical systems are described by 
Hamiltonians of the generic form 



plays the role of g while the hopping amplitude t plays the 
role of r. 

A second example is the constrained XXZ quantum spin- 
1/2 magnet defined on the two-dimensional honeycomb lattice 
by the Hamiltonian-iS 
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where it is assumed that H v can be diagonalized in some pre- 
ferred basis B that spans the Hilbert space TL on which H is 



defined, H commutes with H v , 



Wn^Hjj] — 0, 



H T does not commute with H v 



[H r ,H v ]^0, 



(2.1c) 



(2. Id) 



and the characteristic energy scale U of H v is much larger 
than the characteristic energy scales g and T of H g and H r , 
respectively, 



where the sum over the symbol o is to be understood as a 
summation over all elementary hexagons making up the hon- 
eycomb lattice and the three 2x2 Pauli matrices are denoted 
by <r x , <r y , <r z , respectively. The preferred basis B is the ba- 
sis specified by all the eigenstates with eigenvalues cr| of &\, 
where i runs over all N sites of the honeycomb lattice. In this 
preferred basis the potential energy (I2.3c> and the longitudinal 
one- and two-body interaction J2.3al are diagonal whereas the 
transverse two-body interaction (I2.3b> is not. The exchange 
coupling J, the uniform magnetic field h, and the staggered 
magnetic field h s thus play the role of three different diagonal 
couplings g. 

A third example is the constrained quantum Ising model in 
a transverse field 9,1 ? 



I<?l,|r|«c/. 



(2.1e) 



A famous example of a constrained many-body Hamilto- 
nian is the single-band Hubbard model 7 
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that acts on the fermionic Fock space generated by the cre- 
ation cj CT and annihilation c ia operators for electrons carrying 
the site index i and the spin index a. The preferred basis B is 
the basis specified by all the local fermionic occupation num- 
bers c| CT c io . where i runs over all the lattice sites and a over the 
spin up or down. In this preferred basis, the chemical poten- 
tial ( I2.2ai and the potential energy J2.2d i are diagonal whereas 
the kinetic energy J2.2bl is not. The chemical potential /i thus 



that shares with the spin- 1/2 quantum Hamiltonian (12. 3> the 
same preferred basis. In this preferred basis, the potential en- 
ergy J2.4c> and the one- and two-body interactions (I2.4a> are 
diagonal whereas the one-body transverse field (I2.4bi is not. 
The exchange coupling J, the uniform magnetic field h, and 
the staggered magnetic field h s thus play again the role of 
three different couplings g. 
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Our last example is the quantum Hamiltonian 



For the Hubbard model n = 1 and"i 12 i 13 
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that acts on the bosonic Fock space generated by the bosonic 
creation b\ - +e and annihilation b i i+e operators defined on the 
mid-points of the nearest-neighbor links of the square lattice. 
The preferred basis B is the basis specified by all the local 
bosonic occupation numbers, i.e., the eigenvalues n i i+e of 



acts on the subspace H v of the fermionic Fock space with no 
more than one electron per site. Here, we have introduced the 
fermionic bilinears S i :— c^^^-c^, andfi i := $ ia c ia (sum- 
mation over repeated spin indices is implied). Notice that, at 
half-filling, all contributions to first order in t vanish so that 
r e ff becomes second order in t (n = 2). 

For the constrained XXZ quantum spin- 1/2 magnet on the 
honeycomb lattice ( 12.31 . n = 3 in Eqs. 12.71 whereby the 
constrained Hilbert space is the subspace H Q v defined by all 
states of the form \a\, . . . , a z N ) such that the magnetization of 
each elementary hexagonal plaquette is some integer multiple 
of 3, i.e., 



ft h 
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For the constrained quantum Ising model in a transverse field 
on the honeycomb lattice (I2.4> . n = 6, and the constrained 
Hilbert space is the same as in example d2.3i . 

At last, n = 2 for the bosonic Hamiltonian ( 12.51 whereby 
the constrained Hilbert space is the subspace H v spanned by 



all states of the form 



L (ij) 



.), where (ij) runs over all 



for all the 2N links i, i + e where i runs over the N sites of 
the square lattice and e over its two generating vectors x and 
y. In this preferred basis the potential energy ( I2.5cl and the 
four-body interaction ( I2.5ai are diagonal whereas the kinetic 
energy (I2.5b> is not. The coupling v thus plays the role of g 
while the coupling t plays the role of V. 

By assumption, the contribution H v encodes the largest 
characteristic energy scale U in the problem. As long as all 
other energy scales are much smaller than U, the low-energy 
physics is captured by an effective Hamiltonian H cfi that is 
defined on the Hilbert space restricted to the lowest energy 
eigenstates of Hjj. H cS can be systematically deduced from 
H by treating all the contributions H r to H g r that do not 

commute with H v within perturbation theory. This effective for all N sites z xhis condition is satisfi ed if n 
model is particularly interesting whenever the ground state 
manifold Ha. u of H u is extensively degenerate, as is the case 



2A^ links of the square lattice made of N sites, and the boson 
occupation numbers are restricted to 
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for all examples J2.2l2.5i . 
theory, i? cff is given by 
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and H, 



(n) 



of order zero in T/U, whereby it is understood 

that H cS acts only on the subspace TL Q v of the unconstrained 
Hilbert space with basis B. The order n and the form taken by 
I2.7a> are model dependent. 



H v/u in Ec l 



and only one link out of the four connected to each vertex is 
occupied by a boson. This constrained system is equivalent to 
a square lattice quantum dimer model. 

The last three models will be studied in detail in Sec.lVland 
Sec. lVIl which the reader is referred to for an example-based 
approach. In Sec.llllland Sec. lIVI we now discuss the possible 
topologies of the phase diagram when a strongly constrained 
quantum system possesses a constrained critical regime at in- 
termediary temperatures. 



III. COMPETING CHARACTERISTIC ENERGY SCALES 



From the generic form taken by the low energy Hamilto- 
nian d2.7i acting on the ground state manifold H v and from 
the assumption that |g|, |T| <C U, we deduce the existence of 
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at least three regimes of temperatures provided any one of the 
hierarchies of energy scales 

T«|r cff |,| 5 |«[/, (3.1a) 
|r off |,| 5 | «T«f/, (3.1b) 

|r eff |, \g\ « u « t. (3.ic) 

holds. In addition to these three temperature regimes, strongly 
constrained systems described by the low energy Hamilto- 
nian i2J\ can exhibit a fourth regime if \g\ is much larger 
than all effective off-diagonal couplings. This is the case for 
example when \g\ > \T\ and the off-diagonal terms in H r do 
not contribute to first order in T, i.e., n > 1 in Eq. J2.7bb . 
If so, \T cS | is guaranteed to be much smaller than g and the 
temperature regime < !3.1al can be further subdivided into two 
distinct ones, for a total of four regimes 

T < |r eff | « \g\ < U, (3.2a) 

|r eff | «T« \g\ «(7, (3.2b) 

|r cff |«|. 9 |«T«(7, (3.2c) 

|r eff | < \g\ «£/«T. (3.2d) 

For simplicity, we shall consider only two limiting cases 
that encode the competition between the classical energy scale 
g and the quantum energy scale r eff in this paper. The first 
occurs when \g\ 3> |r eff | in the temperature regime (I3.2c> . 
The second occurs when \g\ <C |r cff | in the temperature 
regime A3. lbl . Both situations can be realized in exam- 
ples ( I2.3I2.5> . In particular, the case of a quantum energy 
scale dominating over the classical one is of relevance to the 
Hubbard model close to half-filling, to quantum XXZ Heisen- 
berg magnets with strong XY exchange anisotropy, to Ising 
magnets subjected to strong transverse magnetic fields, or to 
constrained bosons whose kinetic energy dominates over their 
interactions. We postpone the discussion of the second case 
to Sec. lVllI while we consider here the case \g\ 3> |r eff |. 

The quantum world at T — is governed by a delicate 
competition between the diagonal and off-diagonal energy 
scales. As these characteristic energy scales are varied, quan- 
tum phase transitions can take place between different states 
of matter that usually support long-range order. Universality 
emerges at fine-tuned quantum critical points where the tran- 
sitions between different states of matter are continuous. The 
signature of a quantum critical point can manifest itself as a 
scaling regime as long as temperatures are not too large, say 
when Eq. (I3.2a> holds. 

Upon increasing the temperature from T = 0, one 
leaves the quantum regime ( I3.2ai once T becomes larger 
than the quantum coupling |r off | in the restricted Hilbert 
space. Beyond the quantum regime one distinguishes the three 
regimes d3.2b> . j3.2c> . and ( I3.2dl . If the temperature is much 
smaller than \g\ thermal fluctuations are dominated by the 
classical energy scale g. If the temperature is raised to val- 
ues that are much larger than \g \ but that remain much smaller 
than the characteristic constraint energy U, then the thermal 
fluctuations are predominantly entropic in character, with the 
entropy of the classical constrained phase space isomorphic 



to the preferred basis B QU that spans the ground-state mani- 
fold H u of H v . Finally, once the temperature becomes the 
largest energy scale in the problem, the thermal fluctuations 
are still entropic in character, but now with the entropy of the 
classical unconstrained phase space isomorphic to the basis B 
of H . The transitions between these regimes can take place 
through phase transitions or through crossovers. 

The two entropic regimes ( I3.2cl and (I3.2d> do not always 
need to be qualitatively different. This is the case for the Hub- 
bard model for which all connected spatial correlation func- 
tions between the local electronic densities h ia — $ ia c i(J in 
the entropic regime j3.2ci decay in a qualitatively similar way 
as in the entropic regime l !3.2dt . i.e., exponentially fast with 
separation beyond a characteristic length scale of the order of 
the lattice spacing. The situation may change if one consid- 
ers extended Hubbard models with nearest-neighbor coupling 
V, next nearest-neighbor coupling V', etc, at commensurate 
fillings. We shall comment on this situation in the conclu- 
sions. In examples (I2.3I2.5> . we shall show below that the 
two entropic regimes ( I3.2ci and < !3.2dl are qualitatively differ- 
ent as measured by the temperature dependence and order of 
magnitude of the correlation length characterizing the onset of 
exponential decay in spatial correlation functions. 



IV. THE CONSTRAINED ENTROPIC SCALING REGIME 

In this section we shall study the generic features of the 
constrained entropic regime ( I3.2ci for Hamiltonians of the 
type (12.312.51 . All the assumptions and different cases consid- 
ered here will be supported with explicit examples in Sec.fVl 
and Sec. IVII The reader who may be unfamiliar with the 
physics of constrained models is referred to those two sections 
for an example-based approach. 

Since quantum dynamics is of no qualitative relevance in 
regime (I3.2c> . we shall always assume that |r off |/T — » 
while keeping the ratios \g\/T <C 1 and T/U <C 1 fixed. 
Were it not for the presence of the constraint in regime ( !3.2cK 
set by the large energy scale U 3> T, the condition T ^> \g\ 
would place the system deep into a massive phase, i.e., cor- 
relation functions would decay exponentially in space with a 
characteristic decay length, the correlation length, of the or- 
der of the lattice spacing a. We are going to argue that, in 
examples (I2.3l2.5i . the constraining energy scale U induces 
a correlation length much larger than the lattice spacing and 
possibly increasing with temperature in regime ( !3.2cH 

Without loss of generality g > is assumed in the remain- 
der of this section. 



A. The scaling limit g/T,T/U ^ 

We begin our analysis by considering the scaling limit 
g/T, T/U — > 0. In this limit, all entropic fluctuations are re- 
stricted to the classical configuration space isomorphic to the 
basis B Q v . This hard constraint has dramatic consequences 
on thermal averages in examples J2.3l2.5i . Indeed any spatial 
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spin-spin correlation functions in examples J2.3I2.4I > or spa- 
tial correlation functions between the local bosonic densities 
nuj\ in example (12.51 decay algebraically with separation in 
the scaling regime g/T,T/U — > 0. This is so because ex- 
amples J2.3I2.4I reduce to the non-interacting classical three- 
coloring model that was solved by Baxter in Ref.ll4j, whereas 
example fl2.5i reduces to the non-interacting classical square 
lattice dimer model that was solved by Kasteleyn in Ref . [3 
In either case, it is now understood that the constrained en- 
tropic scaling limit g/T, T/U — > is critical in that correla- 
tion functions decay as power laws in space. 

B. The scaling limit g/T -» 0, < T/U < 1 

In this paper, we shall assume that, if we soften the condi- 
tion that entropic fluctuations satisfy the constraint, i.e., if we 
consider the limit g/T — > holding T/U small but finite, we 
must then impose a cutoff £ ce (T/U) to the algebraic decay in 
space of correlation functions. In effect, we are assuming that 
the operator related to the appearance of defects that violate 
the constraint imposed by the energy scale U is a relevant per- 
turbation to the constrained entropic critical point in the RG 
sense. This assumption is indeed satisfied by examples d2.3i - 
12. 51 but there are no fundamental reasons for it to hold for all 
constrained systems. 

If defects cause the system to flow to a generic uncon- 
strained fixed point with short-range correlations in space, 
the correlation length of the system is controlled by the ratio 
T/U and can be estimated with the following argument. The 
concentration of thermally activated defects that violate the 
constraint is proportional to exp (—a^jU/T) when T <C U, 
where a v > is some non-universal numerical constant spe- 
cific to Hjj, since it represents the fugacity of defects that en- 
ters as the coupling constant driving the system away from the 
constrained entropic scaling limit. Therefore, the correlation 
length of the system is 

4 c (T/L/)~aexp(^-^, (4.1a) 

when T <C U. Here, a is the lattice spacing, d is the dimen- 
sionality of space, and < y v < d is the scaling dimen- 
sion of the operator representing a defect. Upon approaching 
the constrained entropic critical point T/U = 0, the expo- 
nential growth of the correlation length £ CC (T/U) guarantees 
that correlation functions decay in space as power laws over a 
macroscopically large window of length scales 

a«r«4 c (T/[/). (4.1b) 

While the entropic regime ( I3.2cl > of the constrained Hubbard 
model is featureless beyond few lattice spacings, the entropic 
regimes 13.2cl of examples J2.3I2.5I exhibit a scaling behav- 
ior over an exponentially large window of length scales. For 
length scales large compared to the lattice spacing but smaller 
than £, ce (T/U) the system behaves as a non-interacting classi- 
cal system with the constraint fully enforced, whose correla- 
tion functions are captured by a (purely entropic) critical the- 
ory. As the ratio T/U is increased to a value of order 1, the 



correlation length £, CC (T/U) decreases until it becomes of the 
same order of the lattice spacing. 

An interesting possibility occurs if the relevant coupling 
constant that drives the system away from an unstable con- 
strained entropic fixed point generates a correlation length 
£, ce (T/U) that is an algebraic function of T/U as opposed 
to an exponential one in T/U. If so, the correlation length 
£ ce (T/U) ~ (T/U)~ v u for some correlation length exponent 
Vjj. This scaling of the correlation length with temperature is 
identical to that near a quantum critical point, £q CP ~ T -1 / 2 
for some dynamical exponent z. It may thus deceive observers 
expecting quantum critical scaling: the behavior is classical 
as we are in the temperature regime J3.2ci . very far past the 
quantum regime ( 13. 2?\ . 



C. The scaling limit T/U ^ 0,0 < g/T 1 

We now perturb the constrained entropic critical point 
(g/T, T/U) = (0,0) by working at a finite value of g/T 
while T/U = 0. We call scenario I the case when g/T is irrel- 
evant if sufficiently small but becomes relevant beyond some 
finite critical value (g/T) c . We call scenario II the case when 
any finite value of g/T is relevant. Scenario III occurs when 
g/T is exactly marginal up to some critical value (g/T) c , in 
which case the segment < g/T < (g/T) c realizes a line of 
critical points. At last, scenario IV happens when the small 
coupling g/T preserves criticality but strongly alters its na- 
ture; for instance it changes continuously the value of the cen- 
tral charge in the two-dimensional example given in Sec. lVDI 

In the case of scenario I, all properties of the constrained 
entropic critical point (g/T, T/U) — (0,0) survive a suf- 
ficiently small perturbation g/T until it reaches the critical 
value (g/T) c at which a phase transition takes place to a non- 
critical phase of matter selected by the characteristic energy 
scale g. This classical phase transition can be continuous but 
need not be so. The non-critical phase of matter could support 
a conventional classical long-range order such as, say, anti- 
ferromagnetic order. If classical frustration effects are preva- 
lent, the non-critical phase of matter could support less con- 
ventional classical order such as spin-glass order. An exotic 
possibility occurs when classical frustration effects select a 
non-critical phase devoid of any long-range order. 

In the case of scenario II, any finite g/T generates a finite 
correlation length £ u (g/T) by causing the system to order 
into a non-critical phase. This correlation length diverges with 
g/T — > in the case of a continuous classical phase transi- 
tion at infinite temperature, but it may as well remain finite, 
for example if the transition is first order. 

Scenarios III and IV differ from scenario I in that the seg- 
ment < g/T < (g/T) c with T/U = is a line of critical 
points in both cases. The number of critical degrees of free- 
dom is unchanged in case III whereas it does change in case 
IV as a function of < g/T < (g/T) c . 

The constrained quantum XXZ Heisenberg and transverse 
field Ising models ( 12. 3t and ( 12. 4t . respectively, provide ex- 
plicit realizations of scenarios I, II, III, and IV, as we will 
demonstrate in Sec.[V]and Sec. lVIl 
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FIG. 2: (Color Online) - Qualitative phase diagram of the classical constrained system assuming the existence of two phases, a classical 
ordered phase and a disordered one, separated by a phase boundary. The details illustrated in these two figures, such as the phase boundaries, 
only reflect a close neighborhood of the origin of the phase diagram, and not features outside this neighborhood. The dashed lines represent 
curves where g and U are held fixed and only the temperature is varied. These curves originate at the ordered fixed point (oo, 0) at T = and 
end at the disordered and unconstrained entropic fixed point (0, oo) at T = oo. In the limit of g/U <C 1, these curves become infinitesimally 
close to the g/T and T/U semi-axes. The universal physics discussed in this paper occurs in regime I3.2ct . therefore in the region T/U <C 1, 
g/T <C 1 close to the origin of the coordinate system. The shaded region in the Figure represents the scaling entropic region appearing in a 
constrained system due to the proximity to the constrained entropic critical point at the origin. (Left) This diagram encompasses scenarios I, III, 
and IV depending on the behavior of the system along the segment < g/T < (g/T) c with T/U = 0. (Right) Phase diagram corresponding 
to scenario II, where (g/T) c — 0. 



D. Perturbing the constrained entropic critical point with 

< T/U and < g/T 

The fate of the transition at ((g/T) c ,0) on the boundary 
of the phase diagram parametrized by the dimensionless cou- 
plings (g/T, T/U) is model dependent as one moves to the 
interior of the phase diagram. 

For simplicity, we assume that the phase diagram consists 
of two competing phases only. One phase is the basin of at- 
traction of the unconstrained entropic stable fixed point lo- 
cated at (g/T, T/U) = (0, oo). The other phase is the basin 
of attraction of the stable fixed point located at (g/T, T/U) — 
(oo, 0). A cartoon version of this phase diagram is depicted 
in Fig.|5](Left) for scenarios I, III, or IV and in Fig.|5](Right) 
for scenario II. Curves with the dimensionless ratio g/U held 
fixed are represented by dashed lines in Fig. [2] Changing 
the temperature for some given g/U corresponds to moving 
along a dashed line in Fig. |2 The smaller the ratio g/U, 
the larger the temperature range for which the system lingers 
in the vicinity of the constrained entropic critical point (0, 0) 
along a line with g/U held fixed. 

Finally, observe that the location of the phase boundary 
close to the constrained entropic critical point (0, 0) follows 
from 

inia/T) ~ UTIU) (4.2) 

in the case of scenario II, assuming that the phase boundary 
is a line of continuous transitions. One interesting feature of 



scenario II is that the correlation length increases as temper- 
ature increases for a large range of temperatures! The reason 
is that, since the constrained entropic fixed point is unstable 
in both horizontal (g/T) and vertical (T/U) directions, as the 
temperature is raised and one moves in parameter space along 
the dashed line with g/U <C 1, the correlation length first in- 
creases as one approaches the fixed point from the horizontal 
direction, and only starts to decrease as one moves away along 
the vertical direction. The crossover temperature scale is set 
by Eq. ( 14. 2i , and it goes to infinity as U — > oo. 



V. THE CONSTRAINED QUANTUM XXZ HEISENBERG 
AND TRANSVERSE FIELD ISING MODELS 



We are going to illustrate how the high-temperature critical 
scaling picture is realized for the constrained quantum XXZ 
Heisenberg and transverse field Ising models on the honey- 
comb lattice. The effective models are identical for both sys- 
tems, the only difference lying with the order n in T/U needed 

to generate the effective "ring exchange" flpyy m ^1- < !2.7i . 
So it suffices to analyze the constrained transverse field Ising 
model (12. 4> whose effective Hamiltonian ( 12.71 restricted to the 
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subspace H Q v takes the form 



A' 



A 



i=l 



(ij) 



lo«=i ) 



(5.1) 



where the indices i v ,v = 1, . . . , 6 label the sites around an 
elementary hexagon o. The coupling h describes a uniform 
magnetic field, h s describes a staggered magnetic field, and J 
describes an exchange interaction between nearest-neighbor 
sites of the honeycomb lattice. J > favors ferromagnetic 
order while J < favors antiferromagnetic order. Both the 
classical antiferromagnetic and ferromagnetic ground states 
are compatible with the constraint that the magnetization of 
each hexagonal plaquette has to be ±6 or 0. These couplings, 
h, h s , and J, are three g-like couplings that we analyze below. 

The temperature in regime I3.2ct is small compared to 
U and large compared to |r rff | oc |r 6 /Z7 5 |. The approx- 
imation of neglecting either violations of the constraint or 
the off-diagonal part of the quantum Hamiltonian should 
thus be a good starting point. If so, in the scaling regime 
\h\/T, \h s \/T, \g\/T,T/U -> 0, the model reduces to a non- 
interacting Ising model on the honeycomb lattice, with the 
constraint that the magnetization of each hexagonal plaque- 
tte has to be ±6 or 0. This model maps onto the non- 
interacting three-coloring model which was studied by Bax- 
ter and whose entropy can be computed exactly in the ther- 
modynamic limit— He also showed that the model exhibits 
algebraically decaying spatial correlations and as such is crit- 
ical. We shall call this constrained entropic scaling limit 
\h\/T, \h s \/T, \g\/T, T/U -> the Baxter critical point. The 
long-wavelength, low energy limit of this model is captured 
by a conformally invariant field theory with central charge 
c = 2& 

As the temperature is lowered, the effects of the couplings 
J, h, h s and any other couplings compatible with the sym- 
metries of Hamiltonian i5.l\ need to be taken into account. 
That is, we need to decide what perturbations are relevant, 
marginal, and irrelevant at the Baxter critical point. Infinites- 
imally close to it in parameter space one can use perturbative 
renormalization group arguments. At a finite distance away 
these methods fail and one must rely on numerical tools to 
explore the stability of the Baxter critical point. 

Natural choices for perturbations of the Baxter critical point 
are a uniform magnetic field h, a staggered field h s , and a 
nearest-neighbor interaction J. The system is still exactly 
solvable in the presence of the h s coupling alone. In the pres- 
ence of either the h or J coupling, the system is no longer 
exactly solvable, and we chose to resort to numerical trans- 
fer matrix calculations in order to investigate the fate of the 
Baxter critical point. 

Since we are interested in distinguishing between critical- 
ity and any long-range ordered or disordered gapped phase, a 
convenient choice is to measure the central charge of the sys- 
tem. This can be obtained from the coefficient of the largest 
finite-size scaling correction to the free energy of a semi- 



infinite system with periodic boundary conditions in the fi- 
nite directionpiiiS The central charge is known to be strictly 
non-zero if the theory describing the long-wavelength behav- 
ior of the system is critical (massless). Numerically, a mas- 
sive phase is signaled by the vanishing of the measured c, for 
the finite size corrections vanish faster than in a critical sys- 
tem. (There are conformally invariant topological field theo- 
ries with c = 0, but the numerically measured c = is here 
more trivially an indication of a massive phase.) From similar 
calculations on the semi-infinite system one can also obtain 
the scaling dimensions of the operators in the conformal field 
theory describing the long wavelength behavior of the system. 
In addition to providing a better understanding of the critical 
regime, the scaling dimensions are known to either vanish or 
diverge as the system becomes massive, and they can be used 
to confirm the central charge results. 

In order to compute the central charge and scaling dimen- 
sions of the system in a cylindrical geometry, we made use of 
transfer matrix techniques in combination with sparse matrix 
diagonalization routines from the free package ARPACK. Our 
results for the central charge are obtained either as a function 
of h/T or as a function of J/T. We did not consider the case 
of h and J simultaneously present, nor the case when either is 
present together with h s . 



A. Uniform field h - Scenario I 

The numerical results for a uniform magnetic field h are 
presented in Fig. [3] [as the transformation h — > —h and 
a\ — > —of leaves the Hamiltonian unchanged, it is sufficient 
to consider the case h/T > 0]. The dependence on h/T of 
the central charge of a honeycomb lattice wrapped around a 
cylinder is shown in Fig. [3] (Top), for different values of the 
cylinder radius. The dependence on h/T of the two smallest 
scaling dimensions is shown in Fig. [5] (Bottom). 

The numerical calculations proceed in two steps, to be re- 
peated for each value of the reduced magnetic field h / Tjii^ 



First, the three largest eigenvalues A^ '' > Ay > h}-"' of 
the transfer matrix corresponding to the diameters L • = 2j, 
j = 2, . . . , 7 of the cylinder are computed. Here Lj/2 corre- 
sponds to the number of hexagonal plaquettes, or equivalently 
2Lj is the number of spins in a row of the infinite cylinder. 

The largest eigenvalue A^ yields the dimensionless free en- 
ergy per spin 



(2) 



fi 



1 

2L, 



ln(Af) 



(5.2a) 



(1 2) 

while the next two subleading eigenvalues A^ ' ' yield the di 
mensionless excitation energies 



A/f =^.n 



A 



(or 



A 



ik) 



1,2. 



(5.2b) 



Second, a value of the (finite size) central charge c L L is 
obtained from the finite-size scaling fit performed on two con- 
secutive values of the free energy, i.e., on data points of the 
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FIG. 3: (Color Online) - (Top) Behavior of the (finite size) cen- 
tral charge c L L as a function of h/T, obtained from finite-size 
scaling of the free energy computed via transfer matrix. From top 
to bottom, the different curves correspond to increasing diameters 
Lj = 4, • ■ • , 14 for the systems used in the scaling fit to compute 
the value of the central charge c L >L ^. The system in presence of 
a field is symmetric upon the sign change h «-» —h, therefore only 
the h > axis is shown. Notice that the high-temperature criticality 
is robust with respect to a uniform field and it survives at large but 
finite temperatures. (Bottom) Behavior of the two smallest scaling 
dimensions allowed by the conformal field theory as a function of 
h/T. The extrapolated values for L — > oo are shown here for sim- 
plicity. The fact that both (degenerate) scaling dimensions remain 
constant in the critical regime suggests that the whole critical phase 
is described by the same conformal field theory. The lines between 
data points are guides to the eyes. 



tyP e {(Lj,fj), j = i,i + 1}- This is repeated for i = 2 to 
i = 6. The (finite size) scaling dimensions y£ L . k = 1,2, 
follow from finite-size scaling fits on two consecutive values 
of A/j , i.e., on data points of the type {(Lj, A/j^), j = 

What we are after is not so much the finite size values 
c L L and Dj° L as their values in the thermodynamic 
limit L i — > oo (i — > oo). We present finite size data c L L 
in Fig. 13 (Top) to illustrate the finite size corrections while 



we present the extrapolated values j/ fc ) = lirrij—xx, y L L in 
Fig.|3](Bottom). The dependence on h/T of the central charge 



z L L that we deduce from Fig. |3](Top) is 

2, 0< 




< 



(5.3) 



From Fig.[5]one reads that the Baxter critical point is robust 
to the introduction of a uniform magnetic field. According to 
Fig- El (Top), the system remains critical with central charge 
c = 2 over the finite interval < h/T < (h/T) c , before 
entering an ordered phase through a (first-order) phase tran- 
sition at (h/T) c . This is confirmed by the behavior of the 
scaling dimensions shown in Fig.[5](Bottom), which seem to 
rapidly vanish/diverge across the transition at (h/T) c , respec- 
tively. Also, according to Fig. [5] (Bottom), the smallest scal- 
ing dimensions are unchanged along the segment < h/T < 
(h/T) c . These numerical results support the conclusion that 
a small coupling h/T is irrelevant at the Baxter critical point, 
and scenario I is realized along the segment. 

This behavior is perhaps surprising if compared to the effect 
of a uniform magnetic field in the unconstrained Ising model 
on the honeycomb lattice. In that case, the uniform magnetic 
field is a relevant perturbation that causes the system to order 
at any finite temperature. The origin of this difference is due 
to the large depletion of configurations with finite magnetiza- 
tion induced by the projective action of the constraint. The 
entropy of the system as a function of magnetization m seems 
to acquire a cusp at m — that leads to a strong first order 
transition at finite temperature in the presence of a uniform 
magnetic field. 



B. Staggered field h a - Scenario II 



The case of the staggered magnetic field h s has been solved 
exactly by Baxter 14 in the limit T/U — > 0, hjT arbi- 
trary. The model exhibits an infinite order phase transition 
as h s /T — > 0. The staggered field is a marginally relevant 
coupling to the Baxter critical point h s /T, h s /U,T/U — » 
0p& i.e., a staggered magnetic field realizes scenario II of 
Sec. II V (J In particular, this implies that any small staggered 
field h s /T induces an exponentially large correlation length 



£u( h s/ T ) ~ aexp a h 



T 
K 



(5.4) 



where a h is some dimensionless number. 

The divergence of the correlation length ( 15. 41 as h s /T — > 
is cut off at the crossover temperature 



T 



(5.5) 



This estimate follows from the finite correlation length in- 
duced by the fact that T/U, although large, is finite. Indeed, 
the Baxter critical point is unstable to constraint-violating de- 
fects that appear as soon as T/U is finite. The dependence 
on T/U of the corresponding correlation length £, CC (T/U) is 
governed by the most relevant operator that introduces vio- 
lations of the constraint in the spin language. This operator 
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corresponds to the insertion of fractional vortices in the con- 
tinuum theory that describes the Baxter critical point. The 
scaling dimension ify of this operator is given by = 2/9 
according to Ref. 19 in the case when there are no interactions 
added to the model. This is indeed in agreement with our nu- 
merical results for h/T, J/T — > (see Figs. [5] and |4j. The 
estimate (14. lai thus becomes 



L C (T/U) 



aexp I a u 



9U 
16T 



(5.6) 



and the crossover temperature (J53J follows from solving 
Eq. 63 with the help of Eqs. £D and Ol . 



C. Nearest-neighbor interaction J > - Scenario III 

The numerical results for a uniform nearest-neighbor inter- 
action J between the Ising spins are presented in Fig.|4] The 
dependence on J/T of the central charge of a honeycomb lat- 
tice wrapped around a cylinder is shown in Fig. |4] (Top), for 
different values of the cylinder radius. The dependence on 
J/T of the two smallest scaling dimensions is shown in Fig.0] 
(Bottom). 

We consider first the ferromagnetic side of the interaction, 
J > 0. According to Fig.|4](Top), for small but finite J/T, the 
central charge remains larger than or equal to the noninteract- 
ing value c = 2, until the central charge drops to zero abruptly. 
This suggests that the correlation length remains infinite over 
the finite interval < J/T < (J/T)f\ until the system un- 
dergoes a strong first order transition at (J/T)i F ^ . This is con- 
firmed by the behavior of the scaling dimensions, which seem 

( F) 

to rapidly vanish/diverge across the transition at ( J /T) c , re- 
spectively. Moreover, Fig.0](Bottom) suggests that the lowest 
scaling dimensions vary continuously with J/T, even though 
numerics alone cannot be deemed conclusive on this issue. 

The behavior of the central charge in the critical range 
< J/T < (J/T)c F) is different from the uniform field case, 



as c 



is seen to grow significantly upon approaching the 

( F\ 

critical value (J/T) c for a fixed system size This growth 
can be explained by exponential corrections to finite size scal- 
ing that are known to occur at a first-order phase transition. 
Indeed, extrapolating the curves in Fig.|4]in the limit L ■ — > oo 
yields an approximately constant value of c = 2 over the inter- 

( F} 

val < J/T < (J/Tf c The extrapolated curve is however 
rather noisy due to the limited range of numerically accessible 
system sizes and it is not shown here. The reason why these 
corrections are so strong in the case of a nearest-neighbor fer- 
romagnetic perturbation compared to the case of a uniform 
magnetic field deserves further study. 

Our interpretation of Fig. 0] is that, in the thermodynamic 
limit, the central charge is constant on the segment < 

( F) 

J/T < ( J/T)c while the scaling dimensions are not. If so, 

(F) 

the segment < J/T < (J/T) c realizes a line of critical 
points, i.e., J/T is a marginal interaction along this segment, 
thus realizing scenario III. This interpretation agrees with the 
perturbative RG calculation from Ref.|20j 
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FIG. 4: (Color Online) - (Top) Behavior of the (finite size) cen- 
tral charge c L L ^ as a function of J/T, obtained from finite-size 
scaling of the free energy computed via transfer matrix. From top 
to bottom, the different curves correspond to increasing diameters 
Lj = 4, • • • , 14 for the systems used in the scaling fit to compute 
the value of the central charge c L L ^ . Positive values of J corre- 
spond to a ferromagnetic coupling, while negative values correspond 
to an antiferromagnetic coupling. Notice that the high-temperature 
criticality is robust with respect to uniform nearest-neighbor interac- 
tions and it survives at large but finite temperatures both for positive 
and negative J. Antiferromagnetic interactions however deeply af- 
fect the critical behavior, inducing what appears to be a continuously 
varying central charge. (Bottom) Behavior of the two smallest scal- 
ing dimensions allowed by the conformal field theory as a function 
of J/T. The extrapolated values for L — > oo are shown here for 
simplicity. Contrary to the uniform field case, the scaling dimen- 
sions vary on the AF side (J /T) ( C AF) < J/T < 0. Variations of the 
scaling dimensions on the F side < J/T < (J/T)c F ' cannot be 
resolved numerically. The lines between data points are guides to the 
eyes. 



The equilibration properties of the constrained XXZ 
Heisenberg or transverse field Ising model in the ferromagnet- 
ically ordered phase exhibit rather peculiar features, and have 
been discussed in Ref. 9. Based on those results, we expect 
quantum glassiness to appear in the system when the temper- 
ature is lowered across the transition to the F ordered phase, 
at least as long as the transition temperature is small enough 
for the U -violating defects not to play a significant role in the 
equilibration process (Fig.|5|i. 
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D. Nearest-neighbor interaction J < - Scenario IV 



The behavior of the central charge is even more surprising 
on the antiferromagnetic side of the interaction J < 0. The 
values taken by c L ^ in Fig. |4] (Top) are seen to drop be- 
low c = 2 for J/T < —0.1 as soon as the system size is suffi- 
ciently large. As the diameter of the cylinder is increased from 
L i — 4 to L i — 14, the dependence ofc L L on J/T seems 

to indicate the existence of two distinct regimes for J/T. On 
the one hand, the c L L appear to collapse onto a nonvan- 

ishing value of c for sufficiently large J/T (close to J/T = 0) 
in the thermodynamic limit. On the other hand, consecutive 
c L L remain well separated from each other for sufficiently 

small (negative) J/T, thereby suggesting a vanishing limiting 
value c for the central charge in the thermodynamic limit. 



We interpret our finite size simulations as signaling the ex- 
istence of a (continuous) transition at a finite (J/T)i AF ^ ~ 
—0.2, which is in agreement with variational mean-field re- 
sults by Cirillo et al in Ref. |2l[ Our simulations for the two 
smallest scaling dimensions also agree with this interpreta- 
tion. While they are degenerate and they vary continuously 
in the temperature region (J/T)^ F < J/T < 0, they split 
and rapidly vanish/diverge as soon as J/T < (J/T)c , re- 
spectively. As for our numerical results for the central charge, 
we interpret them as indicative of one of three different pos- 
sibilities. The first possibility is that the phase transition at 
(J/T)c F%> separates a phase with vanishing central charge 
below (J/T) AF and a line of critical points with continu- 
ously varying central charge above (J/T)^^ that interpo- 
late between the values c = 3/2 at (J/T)i AF ^ and c = 2 
at infinite temperature. The second possibility is that the fi- 
nite size data when (J/T)i AF " > < J/T < are the signature 
of an infinite sequence of stepwise increases of the central 
charge in the thermodynamic limit that interpolates between 
the values c = 3/2 atf J/T)c AF ^ and c = 2 at infinite tem- 
perature. (See Ref. l22l for a possibly related phenomenon.) 
And finally, the third possibility is that the finite size data for 
c L L collapse in the thermodynamic limit to the central 



charge c = when J/T < (J/T) (AF) and c = 3/2 when 

{J/T)[ AF) < J/T < 0. The plateau with c = 3/2 would 
correspond to a conformal field theory built out of three Ma- 
jorana fermions and endowed with a supersymmetry. More 
accurate simulations, i.e., simulations that can access larger 
diameters of the cylinder on which the honeycomb lattice is 
wrapped, are required to select which of these possibilities 
corresponds to the correct thermodynamic limit. 
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FIG. 5: (Color Online) - Illustration of the phase diagram of the con- 
strained Ising model in the presence of a nearest-neighbor coupling 
J, discussed in Sec.|V| U is held fixed in the figure, and only two out 
of the three parameters of the model, namely J (horizontal axis) and 
T (vertical axis), are shown. The effects of the quantum energy scale 
r become dominant only at very small values of the coupling J, in 
the interval — jF e ff| < J < |r e ff|, and for very small temperatures 
T < |r e ff|, giving rise to the quantum regime shown in the picture. 
At low temperatures, but for \J\ S> |r e ff| the system exhibits two 
classical ordered phases, a ferromagnetic (F) one for J > and an 
antiferromagnetic (AF) one for J < 0. While the transition to the 
AF ordered phase is continuous, the transition to the F ordered phase 
is strongly first order (dotted line). As discussed in Ref. 9, quantum 
glassiness is expected to appear in the system when the temperature 
is lowered across the transition to the F ordered phase, at least as long 
as the transition temperature is small enough for the U -violating de- 
fects not to play a significant role in the equilibration process. 



Possible experimental realization of the constrained Ising 
model in a transverse field 



A possible physical realization of the constrained classical 
Ising model has been discussed in Refs. ll9l20l in the form of 
lattices of superconducting devices with broken time-reversal 
symmetry. Finding an experimental probe mimicking a stag- 
gered magnetic field might allow the observation of a correla- 
tion length that increases with temperature over a large win- 
dow of temperatures in regime d3.2ci . 



VI. THE SQUARE LATTICE DIMER MODEL 



The phase diagram for the constrained Ising model in the 
presence of a nearest-neighbor interaction is summarized in 
Fig.0 



We are now going to illustrate how the high-temperature 
critical scaling picture is realized for the bosonic model in 
example J2.5b . The effective Hamiltonian i2.ll restricted to 
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the subspace H v takes the form 



a3 *-> y) 



i = l 



N 



-r eff ^ J e + i ( a ^+-+ a .+«. 1 +«+-)- i ( a ^») + H.c. 

i=l 

+ ... (6.1) 

Here, a i i+e is the Hermitian operator canonically conju- 
gate to the local bosonic number operator 



i.e.. 



their 

commutator is the C-number i. Under the assumption that 
\v\ ~ \t\ <C U the coupling T e ff oc t 2 /U of the off-diagonal 



term in the effective Hamiltonian i6.l\ is much smaller than v. 
If we relax this assumption by allowing for all possible values 
of the ratio u/r e ff, Eq. J6. It is nothing but the square-lattice 
quantum dimer model introduced by Rokhsar and Kivelson in 
Ref.il 

We are after the high-temperature universal behavior of 
regime (I3.2c> . To this purpose, we consider first the con- 
strained entropic scaling limit \v\/T, T/U — > for which the 
model reduces to the square-lattice classical non-interacting 
dimer model, which was studied by Kasteleyn in Ref.[l5l He 
showed that the entropy can be computed exactly in the ther- 
modynamic limit. He also showed that the model exhibits 
algebraically decaying spatial correlations, and, as such, it is 
critical. We shall call the constrained entropic scaling limit 
\v\/T, T/U — > the Kasteleyn critical point. 

The Kasteleyn critical point is captured by a height model 
which, in the long wave-length limit, is described by the two- 
dimensional conformally invariant field theory 2 ^ 



S = ttK I d 2 x 



with stiffness specified by 



K = - 
2 



and central charge 



(6.2a) 



(6.2b) 



(6.2c) 



A microscopic dimer is represented in the field theory ( 16. 2> 
by a linear combination of two field operators. The first one is 
the charge q e = ±1 "vertex operator" exp (i27rg e 0) where the 
sign assignment has to do with the lattice being bipartite. The 
second one is the "dipolar operator" e • V0 where e is one of 
the two basis vectors of the square lattice such that (i, i + e) 
is the pair of sites covered by the dimer (for some site i). The 
two-point correlation function 



s i27r(<A(x) -<£(»)) 



K 



y\ 



l/K 



(6.3) 



decays with the exponent 2 when K = 1/2. The two-point 
correlation function 



(6.4) 



e • V0(se) e • V0(y) 



\x - y\ 



decays with the exponent 2 for any stiffness K > 0. The 
value ( I6.2b> is thus special in that the scaling dimensions of 
the charge q e = ±1 vertex and dipolar operators are degener- 
ate and equal to 1 . 

A monomer at site i of the square lattice is a defect in the 
dimer covering since site i is not the end point of a dimer. A 
finite concentration of monomers is represented in the field 
theory ( 16. 21 by the local charge q m = ±1 vertex operator 
exp (i2irKq m ip) for the field ip dual to <j>, 



( dip dip \ I dtp d<j> 

\ dx ' dy J \ dy 1 dx 

The two-point correlation function 

,i27ric(^(<B)— fiy)) 



\x - y\ 



K 



(6.5) 



(6.6) 



decays with the exponent 1/2 when K = 1/2. The charge ±1 
monomer is represented by a strongly relevant operator with 
scaling dimension 1 /4. 

In the following we will focus on the temperature region 
where the quantum effects are negligible for the thermody- 
namic properties of the system, i.e., when |r eff | <T< Lf. 



A. The case v = in bipartite lattices - Scenarios II and III 

It is useful to first consider the situation when v = 0, and 
add some other perturbations to dimer models like those re- 
cently studied by Sandvik and Moessner in Ref.|25l Introduce 
a small fugacity for dimers to cover bonds that, although of fi- 
nite length, extend beyond nearest-neighbor bonds. Consider 
separately two cases: 1) when the long bonds break the bipar- 
titeness of the lattice by connecting sites in the same sublat- 
tices, and 2) when the bonds preserve the bipartite nature of 
the lattice. 

We interpret these perturbations within the framework of 
the field theory as follows. The introduction of appropriate 
chemical potentials leads to a small fraction of the longer 
dimers. If one of these longer dimers is simply removed from 
a configuration, one is left with two monomers in the system. 
So a longer dimer can be thought of as two monomers at its 
endpoints that are chained together (see Fig. |6j. In case 1) a 
longer dimer connecting sites in the same sublattice can be in- 
terpreted as a pair of chained monomers with equal charges. 
The operator product expansion (OPE) of the corresponding 
charge q m = ±1 vertex operators leads to charge q m = ±2 
vertex operators. These composite operators have scaling di- 
mension 2 2 x K/2 — 1 when K = 1/2. Thus, they are 
relevant and they open a gap in the system. The underly- 
ing structure consists of a disordered arrangement of dimers 
with short range correlations. Hence, the Kasteleyn critical 
point should be strongly unstable to perturbations that break 
the sublattice symmetry. 26 The simplest perturbation of this 
type is to allow next-nearest neighbor bonds to be covered by 
dimers. Other analytical arguments and numerics are consis- 
tent with this expectation. 25,26 . Therefore, scenario II is real- 
ized upon the introduction of any finite chemical potential for 
second-neighbor dimers in the system. 
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FIG. 6: (Color Online) - Example of longer dimers that can be 
thought of as pairs of monomers chained together. When the end- 
points of a long dimer sit at sites within the same sublattice, the 
chained monomers have equal charges. When the endpoints sit at 
opposite sublattices, the chained monomers have opposite charges. 



In case 2) a longer dimer connecting sites in opposite sub- 
lattices can be interpreted as a pair of chained monomers with 
opposite charges. The OPE of the corresponding vertex op- 
erators (when all the possible directions of the long dimer 
are added) will lead to a |Vr/>| 2 term that renormalizes the 
stiffn ess. Th is is also consistent with the numerical results of 
Refs. l25l26l Scenario III is thus realized. 

We would like to note that violations of the constraint in 
the form of monomers are relevant perturbations that drive 
the system away from the the Kasteleyn critical point, as in 
all the situations discussed is Sec. lIVI The finite correlation 
length j4.^ thereby generated is given by 



477 

t cc (T/U) ~ aexp ( a v — 



(6.7) 



when T -C U. 



B. The case v ^ in bipartite lattices - Scenario III 

We are going to argue on the basis of numerics that scenario 
III can also be realized by perturbing the Kasteleyn critical 
point with the interaction v in Eq. ( 16. U . 

To this end, observe that H v , defined by Eq. line (I2.5ai . 
counts the total number of elementary square plaquettes of 
the square lattice that are flippable in the dimer basis repre- 
sentation. Here, a flippable plaquette is an elementary square 
plaquette that has two occupied edges (dimers). This model 
has already been studied by Alet et al. in Ref.|27 for negative 
values of the coupling constant v. Here we consider both pos- 
itive and negative values, and good agreement with the previ- 
ous results is found where they overlap. 27 



As exact analytical results are no longer available, we use a 
numerical approach similar to the one described in Sec.|V] We 
compute the central charge of the system as well as the scal- 
ing dimensions of two specific operators via transfer matrix 
techniques, using finite size scaling fits. The accessible sys- 
tem sizes are Lj = 2j, j = 3, . . . , 8 where Lj is the number 
of square plaquettes across the periodic direction of the sys- 
tem, or equivalently 2L, is the number of edges. As discussed 
by Alet et al in Ref. 27| the critical regime of this system is 
captured by a c = 1 two-dimensional conformal field theory 
of the Coulomb gas type, with continuously varying stiffness. 
They also computed the scaling dimensions of some known 
operators in the conformal field theory, namely those corre- 
sponding to the electric and magnetic vortices in the Coulomb 
gas picture. From these measurements one can then obtain the 
values of all other scaling dimensions present in the conformal 
field theory and discuss the nature of the critical phase and of 
the phase transitions. 

The central charge as a function of the coupling constant 
v/T is plotted in Fig. (Top). Thanks to the good conver- 
gence in the L — > oo limit, we report the extrapolated values 
of c instead of each separate curve for increasing system size. 
From the existence of a c = 1 plateau 



^jij (columnar) < „/ T < (v/T) ( ^ S: 



gered) 



(6.8) 



that extends to the left and right of the Kasteleyn fixed point 
v/T,T/U — > 0, we infer that criticality is preserved upon 
introducing a small but finite coupling v/T. Moreover, the 
behavior of the scaling dimensions presented in Fig. (Bot- 
tom) tells us that the stiffness K in the critical theory ( I6.2al > 
varies continuously along the line of critical points. Therefore 
we conclude that the introduction of the coupling v realizes 
scenario III. 

The transition at T = j^ column<ir ) nas been characterized by 
Alet et al. The quasi-long-range ordered phase ( 16.81 under- 
goes a Kosterliz-Thouless transition at («/T)c ColllIimar to an 
ordered phase. This phase is characterized by an alignment of 
parallel dimers along columns or rows and is therefore called 
the columnar phase. 

On the other side, for v > 0, our results show that the quasi- 
long-range ordered phase d6.8l > terminates at T — j , c ( sta ss ered ) 
where it undergoes a first-order phase transition to an ordered 
phase. This phase is characterized by a staggering of parallel 
dimers along two consecutive columns or rows and is there- 
fore called the staggered phase. 

The phase diagram in the regime ( I3.2cl has the same topol- 
ogy as the one for the constrained Ising model in the presence 
of a nearest-neighbor interaction (Fig.[5j. 



VII. QUANTUM DOMINATED REGIME, \g | < |r efI | 

In Secs. lIIIHVII we assumed that |r eff | was the smallest en- 
ergy scale in the problem. This is to be expected whenever 
r cfr = T(T /U) n ~ 1 is highly suppressed because of the order 
n needed for virtual processes to return to an allowed config- 
uration. Even if we initially had \g\ < \T\, it is likely that 
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FIG. 7: (Color Online) - (Top) Behavior of the extrapolated central 
charge as a function of v/T, obtained from finite-size scaling of the 
free energy computed via transfer matrix. (Bottom) Behavior of the 
scaling dimensions corresponding to the electric and magnetic vertex 
operators as a function of v/T, also obtained from finite-size scaling 
arguments and transfer matrix calculations. Observe that the scaling 
dimensions of the electric and magnetic vertex operators are 1 and 
1/4, respectively, when v/T = 0. The lines between data points are 
guides to the eyes. 



r off | < |g| < |r|, since |T| <C U. However if the coupling 
constant \g\ for those terms that commute with the constraint 
is small compared to |r cff |, or in particular, if g = 0, the dis- 
cussion above must be revisted. 

Let us consider here the case g = for simplicity. Indepen- 
dently of whether r eff is suppressed with respect to T or not, 
there are now only three regimes of temperatures: 



t« |r eff | «[/, 
|r cff | «r«[/, 
|r cff |«f/«T. 



(7.1a) 
(7.1b) 
(7.1c) 



The constraint, as well as all other energy scales in the sys- 
tem, becomes negligible in regime ( 17. lcl and the resulting 
physics is that of a featureless high-temperature phase, say a 
paramagnetic phase if the degrees of freedom are exclusively 
magnetic. 

As the temperature is lowered down to regime ( 17. lbi . the 
system is still classical (at least from the point of view of its 
thermodynamic properties) but the constraint is now enforced. 



The physics is controlled by the proximity to the constrained 
entropic critical point at T cff |/T, T/U — * 0. 

Finally, in regime 17 . 1 ai the system is fully quantum me- 
chanical. 

When \g\ is finite but much smaller than \T cS \, these argu- 
ments should still be valid. Therefore, we can draw a qualita- 
tive phase diagram for a generic strongly constrained quantum 
system at fixed g as a function of r off . In Fig. [8] we represent 
the phase diagram of a system exhibiting a zero temperature 
phase transition at a finite value of the ratio T cff /g. Notice 



T~lg| 




FIG. 8: (Color Online) - Generic phase diagram for a strongly con- 
strained quantum system that satisfies conditions 1) to 6) from Sec.0 
Parameter space encodes the competition between two energy scales, 
the temperature T and a characteristic energy scale r cff given by the 
effective coupling of the off-diagonal term in the Hamiltonian, while 
g and U are held fixed. Here we chose to represent the case of a sys- 
tem exhibiting a zero temperature phase transition at a finite value 
of the ratio T cS /g. A quantum critical scaling regime, if it exists, 
is restricted to a region represented by the upper half of a disk of 
radius r cff and centered at the origin (g, T) — (0, 0) of parameter 
space. As in the case of Fig. Q the system exhibits a constrained en- 
tropic scaling regime at temperatures intermediate between the small 
characteristic quantum energy scale r cff and the large characteristic 
energy scale U set by a strong constraint. At a fixed temperature, the 
constrained entropic scaling regime terminates in a phase transition 
that needs not be continuous upon increasing |r eff |. At fixed F eff , the 
constrained entropic scaling regime crosses over to the conventional 
high-temperature phase, say a paramagnetic one for spin degrees of 
freedom, when T is of the order of U. 



that the shaded dome around the origin or parameter space 
contains both the classical ordered phase that onsets when 
\g\ becomes larger than |T eff |, and the quantum critical scal- 
ing regime that is expected to appear in a cone-shaped region 
above the zero temperature quantum critical point. The actual 
details of this shaded region are highly system specific. 

It is worth mentioning that the existence of a zero tempera- 
ture phase transition at r eff ~ ±\g\ is by no means a necessity. 
It is also possible - for example - that the coupling T cff never 
destroys the ordered phase determined by the diagonal cou- 
pling g at zero temperature, e.g., if it favors the same type of 
order. In this case there would be no quantum gapped phase 
in Fig. [8] and the dome above the origin of parameter space 
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would stretch out from (—00, 0) to (00, 0), without any phase 
transition at (±00, 0). The phase under the dome would be 
uniform and determined by the (fixed) value of g. This is so, 
for example, in the quantum Hamiltonian (16. 1> for any fixed 
value of g < 0, where the columnar ordered phase is stable at 
all values of T eff according to the numerical results in Ref.Esl 

VIII. CONCLUSIONS 

In this paper we have presented a mechanism that leads to 
a large temperature regime where critical scaling behavior ap- 
pears as a consequence of having a large energy scale U in 
the problem. Hard constraints are imposed on the system in 
the limit U — > 00, that project the original Hilbert space onto 
a space of allowed configurations satisfying the constraints. 
Critical behavior in the U — > 00 limit occurs when corre- 
lations functions, calculated as uniform averages over all of 
these projected states, are algebraically decaying in space. We 
call the limit U — > 00 followed by T — > 00 at which critical- 
ly emerges a constrained entropic critical point. For large but 
finite U, and in the presence of other couplings \g\, \T\ <C U 
in the problem, the physics at temperatures large compared 
to all other couplings but small compared to U is still con- 
trolled by the proximity to the constrained entropic critical 
point. Indeed, this hierarchy among the couplings opens a 
rather large window of temperatures for which a constrained 
entropic critical regime exists as a consequence of the prox- 
imity to the constrained entropic critical point. We expect the 
constrained entropic critical regime to be qualitatively differ- 
ent in general from the quantum critical regime associated to 
a putative quantum critical point at zero temperature. The pic- 
ture is summarized in Fig. [2 

The mechanism for criticality at high temperature discussed 
here should be applicable to some, but not all, problems with 
a large dominant energy scale. Examples can be found in 
frustrated magnets, say pyrochlore antiferromagnets where a 
magnetic field that induces a magnetization plateau at half the 
value of full magnetization supplies a strong constraint. 29 - 30 ' 31 
Furthermore, since the order n of a virtual process leading to 
the effective kinetic quantum ter m T „ = T(T /U) n ^ 1 is usu- 
ally large (it is of order 9 in Ref. l30l) . the window of temper- 
atures for which the system is in the quantum regime can be 
extremely small (T < |r off <C |T|). Hence, in practice, even 
for the smallest temperatures accessible experimentally, these 
constrained systems should either order or be within the con- 
strained entropic critical regime. What determines whether 
the system orders or not is the presence of another energy 
scale g, set by the coupling strengths of additional terms in 
the Hamiltonian that commute with the U -term imposing the 
constraint. 

Another example for which the idea of high-temperature 
criticality may apply is that of the fluctuations about a so- 
called d-density wave (DDW). 32 In this instance, the order 
may form locally at some large energy A d , but the current 



loop directions may fluctuate. 6 ' 33 At low temperatures, the 
system should order in a given current pattern. Above the 
global ordering temperature, this system should display high- 
energy constrained entropic criticality, because it resembles 
an ice model as long as T is below A d . 

In the single-band Hubbard model, with a strong local on 
site repulsion only, there is no constrained entropic critical 
point in the U — > 00 limit. It has been recently proposed 
by Philip Phillips and co-workers that the order in which the 
limit U — ► 00 and the thermodynamic limit L — > 00 (L being 
the linear size of the system) are taken has implications for 
hole transport in the single-band Hubbard model. 34 The issue 
is that, for finite U, there is a characteristic length scale asso- 
ciated with the distance between doubly occupied sites, and 
hence for system sizes greater than this distance one has a fi- 
nite density of such "defects". The appearance of this extra 
length scale was suggested as a way to resolve the issue of 
the breakdown of the one-parameter scaling picture for quan- 
tum criticality in the cuprates. 35 While we do not know how 
to connect their results to ours, it seems that there is one com- 
mon theme, that high-energy terms find their way to affect the 
physics at intermediate temperatures. 

The chances to find a constrained entropic critical point 
improve if one adds nearest-neighbor, next-nearest neighbor, 
etc., couplings, i.e., in extended versions of the Hubbard 
model, at some commensurate fillings. For example, the pres- 
ence of very strong nearest-neighbor interactions at 1/4 filling 
leads to a classical checkerboard configuration, where one of 
the two sublattices of the bipartite square lattice is filled and 
the other empty. If one now starts changing the doping away 
from 1/4 filling, extra holes or particles will tend to cluster 
in s tripesi2& 3 i2&22i22i 4 iii 2 i Other longer range couplings that 
commute with the constraint could lead to charge ordering, 
perhaps in the form of static stripes or other patterns at low 
temperatures. Now, at high temperatures, when charge order 
is destroyed, these stripes meander and fluctuate, and the sys- 
tem could be in an entropic critical regime. Certainly, these 
are speculative thoughts, in contrast to the concrete cases of 
the constrained quantum Heisenberg and transverse field Ising 
models that we presented as examples of constrained entropic 
criticality. But the discussion suggests how it is not implausi- 
ble that a constrained entropic critical point may play a role in 
the physics of the extended Hubbard model near commensu- 
rate fillings. 
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